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NATURAL VIBRATION AND FLUTTER OF ELASTICALLY SUPPORTED 
CORRUGATION-STIFFENED PANELS - EXPERIMENT AND THEORY 

By Walter L. Heard, Jr., and Herman L. Bohon 
Langley Research Center 

SUMMARY 

The natural vibration and flutter characteristics of two corrugation-stiffened panels 
were studied experimentally and theoretically to determine the effects of finite deflec- 
tional, rotational, and torsional stiffnesses of the boundary supports. Flutter tests were 
conducted at a Mach number of 3 in the Langley 9- by 6-foot thermal structures tunnel. 

The test specimens were restrained against in-plane thermal expansion in the streamwise 
direction; thus, uniform in-plane loading was induced by aerodynamic heating. The exper- 
imental flutter data are compared with flutter theory for orthotropic panels (developed in 
the present report) which includes the effects of finite support stiffnesses as well as 
uniform in-plane loading of the panel. The panel support stiffnesses (determined from 
experimental natural frequencies) and an assumed state of stress are used in the theory 
for comparison with the flutter data. The investigation indicates that finite deflectional, 
rotational, and torsional stiffnesses of the boundary supports must be accounted for in the 
theory in order to match experimental natural frequencies of the test panels. However, 
it is concluded that structural damping (neglected in the theory) strongly influenced the 
flutter behavior of the test panels. Furthermore, the theory shows that the addition of 
torsional stiffening at the supports is a definite asset for increasing flutter dynamic pres- 
sure and may be extremely important in the design of panels for application at hypersonic 
speeds where free thermal expansion may be desired. 

INTRODUCTION 

Current designs of hypersonic vehicles indicate that exterior surfaces may consist 
essentially of corrugation-stiffened panels with flexible boundary supports. Such supports 
may be designed to allow essentially free in-plane thermal expansion of the panels; thus, 
buildup of large thermal stresses due to boundary restraints can be prevented. Panel 
flutter characteristics, however, may be adversely affected by flexible supports unless 
careful attention is given to the design of the supports. (See ref. 1.) 

In reference 2 early experimental flutter data for corrugation-stiffened panels are 
compared with theory for orthotropic panels; the experimental values for flutter dynamic 
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pressure are found to be very much lower than those predicted by the theory. This dis- 
crepancy is attributed to the assumption, in the theory, of infinite deflectional stiffness of 
the supports at the ends of the corrugations. Reference 3 shows that if the supports are 
assumed to have finite deflectional stiffness, the agreement of theory with experiment is 
greatly improved; nevertheless, significant differences remain (ref. 4). In the theory of 
reference 3, however, torsional support stiffness is assumed to be zero, rotational sup- 
port stiffness is assumed to be either zero or infinite, and panel in-plane loading is 
neglected. Thus, the objectives of the present paper are (1) to provide experimental 
natural vibration and flutter data when the effects of the boundary support stiffnesses and 
in-plane loading of the panel may be important considerations and (2) to provide a more 
complete orthotropic -panel flutter theory by adding effects of finite torsional and rota- 
tional stiffnesses and uniform in-plane loading to the theory of reference 3. 

Experimental flutter characteristics are presented for two thermally stressed 
corrugation-stiffened panels. Both panels were tested with the corrugations alined normal 
to the airstream. The leading and trailing edges of the panels were essentially simply 
supported. The side (streamwise) edges of the panels were attached to finite -stiffness 
supports designed to permit simple calculations of deflectional and rotational stiffnesses. 
The torsional stiffnesses of the supports were determined by correlating measured natural 
frequencies of the panels with theory. After the natural frequencies were measured, 
flutter tests were conducted at a Mach number of 3 in the Langley 9- by 6-foot thermal 
structures tunnel. The test specimens were restrained against thermal expansion in the 
streamwise direction; thus, uniform in-plane loading was induced by aerodynamic heating. 
Results are presented in the form of plots of the critical dynamic -pressure parameter as 
a function of change in panel temperature. 

The flutter theory for stressed orthotropic panels presented herein accounts for the 
finite deflectional, rotational, and torsional stiffnesses of the boundary supports at the 
ends of the corrugations. The other two edges are assumed to be simply supported. The 
effects of the various spring supports are discussed and theoretical flutter mode shapes 
are presented. The experimental data are compared with flutter boundaries predicted by 
the theory by using an assumed state of stress. 

SYMBOLS 

The units used for the physical quantities in this paper are given both in the U.S. 
Customary Units and in the International System of Units (SI). Factors relating the two 
systems are given in reference 5, and those used in the present investigation are pre- 
sented in appendix A. 
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B y,m 
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CmjDm 


amplitude coefficient 
quantity defined in equation (B9) 
panel length (x-direction) 
frequency parameter (see eq. (B8)) 
panel width (y-direction) 
constants of integration 
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panel bending stiffnesses in x- and y-directions, respectively 
panel twisting stiffness 
width of beam spring 
Young's modulus 


D 


xy 
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E = EI 


hD2 


f 

h 

I 


frequency 

distance between beam springs (see fig. 2(a)) 
moment of inertia of spring cross section 


K Dj Kr,K T deflectional, rotational, and torsional spring constants, respectively, per unit 
length 

K d ,K r ,K t nondimensional deflectional, rotational, and torsional spring constants, 
respectively (see eqs. (B13)) 
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N x ,Ny 

Nx,cr 
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AT 
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Cr 
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effective length of beam spring 
Mach number 

number of half-waves in streamwise direction 

uniform in-plane loads per unit length in x- and y- directions, respectively 

in-plane buckling load at flutter 

dynamic pressure of airstream 

dynamic pressure of airstream at flutter 

change in average temperature of panel cover sheet 

thickness of support spring 

lateral deflection of panel 

Cartesian coordinates of panel 

function describing shape of natural mode of vibration in y-direction 
(see eqs. (B6) and (BIO)) 

mass per unit area of panel 

parameters associated with characteristic roots (see eqs. (Bll)) 

dynamic -pressure parameter, ■ ■ 

Di^m2 - 1 

dynamic-pressure parameter at flutter 
Poisson's ratio 

Poisson's ratio associated with curvature in y- and x- directions, respectively 
time 
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circular frequency 


io 

co r fundamental frequency of simply supported beam (radians per sec), 

(7r2/a 2 )^Djyy 


Subscript: 

ref reference value 


TEST APPARATUS 
Panels 

Two corrugation-stiffened panels were constructed for the present tests according 
to the construction details shown in figure 1. The panels were identical except that for 
panel I, a/b = 1/2, and for panel n, a/b = 1. The width b of each panel was 20 inches 
(50.8 cm). The panels, consisting of a corrugated sheet and a flat cover sheet (external 
skin), were fabricated from Rene 41. The corrugated sheet thickness was 0.010 inch 
(0.025 cm). The flat cover sheet 0.007 inch (0.018 cm) thick, was spotwelded at approx- 
imately 0.25-inch (0.64-cm) intervals along the length of each corrugation (two rows of 
welds per corrugation). The cover sheet provided a relatively smooth aerodynamic sur- 
face. A Rene 41 doubler strap, 0.60 inch (1.52 cm) wide and 0.030 inch (0.076 cm) thick 
(consisting of three 0.010-inch-thick (0.025-cm) straps), was seamwelded to each of the 
streamwise edges of the corrugation sheet prior to forming; the sheet and straps were 
then formed simultaneously. The straps, which were treated as part of the boundary 
supports and not as part of the panel, strengthened the spring-supported edges of the 
panel for mounting purposes and provided a source of torsional stiffness. At the leading 
and trailing edges, formed aluminum channel section supports were attached to the panels 
as shown in detail A of figure 1. These channels were 1.70 inches (4.32 cm) deep and 
0.050 inch (0.127 cm) thick. Each channel was attached to the panel by a single row of 
rivets. The panel stiffness parameters D^, Dg, and Dxy (identical for the two panels) 
were calculated by use of the method in reference 6 and are given in table I. 

Support Springs and Mounting Arrangements 

The support springs and the mounting arrangements for the streamwise edges and 
the leading and trailing edges are shown in figure 2. Two sets of springs (fig. 2(a)) were 
machined from blued shim stock so that when the panels were mounted, the ends of each 
corrugation of the panel rested on a small cantilever beam spring. The spring thickness 
was 0.015 inch (0.038 cm) for the rectangular panel, and 0.029 inch (0.074 cm) for the 
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square panel. The width of each beam spring was 0.50 inch (1.27 cm); the effective 
length l was assumed to be 1.26 inches (3.20 cm), the distance between center lines of 
the attachment holes. The spring spacing h was 1.00 inch (2.54 cm). A view of a 
streamwise edge of a panel showing the corrugations, doubler strap, and hole for attaching 
the spring is shown in figure 2(b). Figure 2(c) shows the streamwise-edge support 
arrangement. Small aluminum blocks were used as spacers between the spring and the 
panel corrugations to aline the flat surface of the panel flush with the aerodynamic surface 
of the test-section panel holder. These blocks were chamfered as shown to minimize 
changes in rotational restraint of the springs between upward and downward movement of 
the panel edges. Each beam spring was attached to its corresponding corrugation by a 
single machine screw and an aluminum block which was drilled and tapped and inserted 
into the ends of the corrugations. The support arrangement for the leading and trailing 
edges is shown in figure 2(d). 

The panels were mounted in a flat-sided steel panel holder (figs. 3 and 4), which 
extended vertically through the test section. The panel holder has a beveled half -wedge 
leading edge with a cavity on the nonbeveled side 5.0 inches (12.7 cm) deep, 29.00 inches 
(73.6 cm) in the direction of airflow, and 30 inches (76 cm) high. The cavity, which is 
located 26.50 inches (67.3 cm) downstream from the leading edge, is used for accommo- 
dating test specimens. The panel holder is equipped with pneumatic, vertically operating 
sliding doors which cover the cavity area to protect the test specimens during tunnel start 
and shutdown. Aerodynamic fences attached to the doors insure essentially free-stream 
flow conditions over the cavity area. The pressure inside the cavity and behind the test 
specimen is controlled by a vent-door arrangement on the side of the panel holder opposite 
the panel. The test panels were attached to a mounting frame which was inserted in the 
cavity. A test panel mounted in the panel holder in the wind-tunnel test section is shown 
in figure 4. The view is from upstream and the protective doors are shown fully opened. 

Instrumentation 

Five iron-constantan thermocouples spotwelded to each of the panel cover sheets 
at the locations shown in figure 5 were used to measure cover-sheet temperatures during 
the flutter tests. Three variable reluctance-type deflectometers were located behind the 
panel in the panel-holder cavity to detect panel motion and to obtain flutter frequencies. 

The locations of the deflectometers are also indicated in figure 5. All wind-tunnel tests 
were monitored with high-speed 16-mm motion pictures to detect panel motion. The 
exposed surface of each panel was painted with a grid for photographic purposes. (See 
fig. 4.) 

Static pressures in the tunnel and at several locations on the panel holder and in the 
cavity behind the panel were measured by quick- response strain-gage-type pressure 
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transducers. Tunnel stagnation pressures were obtained from static pressures measured 
in the settling chamber. Stagnation temperatures were measured by total-temperature 
probes in the test section. All pressure and temperature data were recorded on magnetic 
tape. Deflectometer output, which was monitored during each test, was recorded on 
oscillographs. 


Wind Tunnel 

The tests were made at a Mach number of 3 in the Langley 9- by 6-foot thermal 
structures tunnel, an intermittent blowdown wind tunnel which exhausts to the atmosphere. 
A heat exchanger is preheated to provide stagnation temperatures ranging from approx- 
imately 300° to 660° F (422 to 622 K) and dynamic pressures ranging from 1500 to 
5000 lbf/ft 2 (72 to 239 kN/m 2 ). Run times are limited to less than a minute. For more 
details regarding this facility and its operation, see reference 7. 

TEST PROCEDURE 

Natural frequencies of each panel were measured immediately prior to each flutter 
test. All natural-vibration tests were conducted at ambient temperature and pressure 
in the tunnel test section with the panels mounted in the panel holder. The panels were 
vibrated with an air jet shaker, which is described in detail in reference 8. 

All flutter tests were conducted at a Mach number of 3.0, dynamic pressures 
ranging from 1500 to 5000 lbf/ft 2 (72 to 239 kN/m 2 ), and stagnation temperatures ranging 
from 303° to 452° F (424 to 506 K). For each test, the stagnation temperature was main- 
tained essentially constant; whereas the dynamic pressure was often varied in an attempt 
to either start or stop flutter. The usual test procedure consisted of establishing test 
conditions and then exposing the panel to the airflow by opening the protective doors. The 
initial test conditions were maintained for an established period of time (usually about 15 
to 20 seconds) to allow the thermal stresses resulting from aerodynamic heating of the 
panel to initiate flutter. The dynamic pressure was then lowered in an attempt to stop 
flutter. If, however, flutter did not occur, the dynamic pressure was increased until 
flutter was observed and then lowered again (time permitting) to obtain a flutter stop point. 
The maximum duration of a test was approximately 45 seconds. 

RESULTS AND DISCUSSION 
Experiments 

Vibration . - Immediately prior to each flutter test, the natural frequencies corre- 
sponding to 1 half-wave in the cross-stream direction and m half-waves in the 
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streamwise direction were measured. The measured frequencies for each mode were 
averaged for all the tests of a given panel, and the average values and deviations from the 
average are given in table II. 

Flutter .- Four flutter tests of each panel were made for this investigation. During 
each test the pressure inside the panel-holder cavity was controlled so that the differ- 
ential pressure across the panel was normally very close to zero and never more than 
±0.01 lbf/in2 (69 N/m2). Flutter was observed in all tests except one test for panel I- 
Pertinent tunnel and panel data at flutter for all tests are listed in table HI. The data are 
also presented in figure 6 in plots of the nondimensional flutter parameter A cr as a 
function of the change in average temperature of the panel. The open symbols in figure 6 
represent the start of flutter; the solid symbols represent flutter stop points, which were 
obtained by decreasing dynamic pressure after flutter was established. For a typical test, 
a tunnel stagnation pressure was preselected. The tunnel was started and was allowed to 
reach steady flow conditions. The protective doors of the panel holder were then opened, 
and the panel was exposed to the flow and aerodynamic heating. Thermal expansion of the 
panel was restricted by the boundary supports; thus as the panel temperature increased, 
thermal stresses were built up and eventually caused the panel to flutter. The dynamic 
pressure was then decreased (which, for two tests of panel I, was also accompanied by a 
slight decrease in panel temperature) until flutter ceased or the test was terminated. 
Flutter stop points were not obtained in all tests, however, because of the time limit 
involved for the blowdown wind tunnel. 

Segments of the flutter boundary for each panel are represented by the solid curves 
faired through the most conservative experimental flutter start points shown in figure 6. 
For comparison with theory, the flutter dynamic pressure at zero in-plane load is deter- 
mined by a linear extrapolation of the curve to AT = 0. This extrapolation procedure 
has often been used with good results for aerodynamically heated isotropic and orthotropic 
panels. As can be seen in the figure, the linear extrapolation yields A C r = 4000 for 
panel I and Acr = 20 000 for panel n. However, using the linear extrapolation procedure 
on the data obtained in the investigation probably does not give the correct results, as will 
be discussed subsequently. 


Comparison With Theory 

In order to have a more complete theory for comparison with the experimental data, 
the flutter theory of reference 3 for orthotropic panels (simply supported at the leading 
and trailing edges) is modified to include a torsional spring constant lOp and a rotational 
spring constant along the lateral edges, and uniform in-plane loads N x and Ny. 

The panel configuration and the spring support system are shown in figure 7; the details of 
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modifications to the theory are presented in appendix B. Converged vibration and flutter 
results were obtained by utilizing up to 40 modes in the analysis. 

Spring constants .- Calculated nondimensional deflectional and rotational spring 
constants Kj) and Kr (see eqs. (B13)), respectively, were obtained by simple strength- 
of-materials methods applied to the idealized structure shown in figure 8. Note that the 
eccentricity of the beam springs to the panel neutral surface was neglected. The elements 
representing the beam springs had a stiffness corresponding to El and were assumed to 
be clamped at the boundaries. The element representing the panel had a stiffness corre- 
sponding to D2 per unit length. By applying conditions of continuity of deflection, slope, 
moment, and shear at the spring-panel juncture, the following expressions for Kp> and 
K R were obtained: 


77 12E /V 3 

K D = ~ T 

7T J 


3£ + l\ 

- b — +11 


— ,-+ II 


L \E b 


^ E b 


Kr = 


7 r l 
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( 2 ) 


where E = 

hD2 

By using these equations, the results are = 0.918 and K R = 0.098 for panel I 
and K D = 4.33 and K R = 0.081 for panel II. 

Because of the nature of the support system, the nondimensional torsional spring 
constants are difficult to evaluate. Therefore, these constants were obtained from the 
measured natural frequencies. However, a lower limit of K-p is obtained by taking into 
account the torsional restraint of the beam spring only. Thus, 


: = Et3bd 

T 6(1 + ju)hD 2 ^ 


(3) 


The lower-limit values are Kx = 0.005 


for panel I and Kx = 0.037 


for panel H. 
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The nondimensional theoretical natural frequencies are compared with experimen- 
tally obtained results in figure 9. The dashed curve in figure 9(a) is the theoretical fre- 
quency based on the calculated values of K R and K R given in equations (1) and (2) and 
the lower limit of Kp. As can be seen, the theory does not agree with experiment. 
However, as shown by the solid curve, when Kp is arbitrarily increased to a value 
of 1.1, theory is brought into excellent agreement with experiment. Since a reliable 
calculated estimate of K T is not available, the higher value of K T shown in the figure 
is presumed to be the more realistic value and should be used in the flutter calculations. 
Thus the final support spring constants for panel I are = 0.918, K R = 0.098, and 

k t = 1.1. 

The spring constants for panel n were determined with the aid of figure 9(b). As 
shown by the dashed curve, the use of the lower limit of Kp again resulted in poor 
agreement between theory and experiment. However, unlike panel I, no value of Kp 
could be found to produce satisfactory agreement of theory with experiment when the 
original calculated estimates of K D and K R were used, as shown by the dot -dash curve 
which corresponds to Kp = °°. Thus, in order to match the experimentally obtained 
natural frequencies, other values of K R and/or K R were required. An expression for 
determining new values of K R and K R was derived by considering the beam spring 
only. For a cantilever beam clamped at the root and elastically restrained and loaded at 
the tip, the deflectional stiffness is related to the rotational stiffness according to 
strength-of-materials theory by 


_ 12E/b\ 3 

D Till 


/. 


+ 1 


k r I + 1 

—71 + 1 

E b 


(4) 


where 0 5 Kr i °°. 

The actual amount of rotation at the beam tip was not measured. However, by 
selecting values of K R , corresponding values of K D were determined from equation (4). 
Then, for each combination of K D and K R , various values of Kp were tried until the 
theory gave the best agreement with the measured natural frequencies. The best agree- 
ment for panel II (k d = 7.1, K R = 2.5, Kp = 8^ is shown by the solid curve in figure 9(b). 
The final values of Kp for panels I and II are seen to be over two orders of magnitude 
higher than those values obtained from equation (3) for the beam only. Thus it would 
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appear that the primary torsional stiffness K-p results from the doubler strap. The 
selection of the final values of K^, Kj^, and K,p were based on the best agreement of 
theory with experiment; these values are given in table IV. 

Flutter predictions for zero in-plane loading .- The selected (final) values of Kp>, 

Kr, and K-p given in table IV were used to define the natural frequencies and corre- 
sponding natural mode shapes for the flutter analysis of appendix B. The variation of 
panel frequencies with dynamic-pressure parameter A is shown in figures 10(a) and 
10(b) for panels I and n, respectively. The behavior of only the first four modes is shown. 
As the dynamic pressure increases, modes two and three approach each other until they 
coalesce. At this point (the maximum point on the loop) the panel becomes unstable and 
flutter occurs for any increase in A. The resulting flutter mode shapes, then, are defined 
predominantly by modes two and three. These theoretical flutter mode shapes are shown 
in figures 11(a) and 11(b) for panels I and n, respectively. 

The critical flutter dynamic-pressure parameter predicted by theory for zero 
in-plane loading is A cr = 11 235 for panel I and A cr = 59 700 for panel n. These values 
of A cr are nearly a factor of 3 greater than the values obtained from the linear extrap- 
olation of the experimental data in figure 6. 

In light of the apparent discrepancy between theory and the extrapolated experimental 
flutter dynamic pressure q cr of each test panel, the question arises about the sensitivity 
of q cr to variations in panel support stiffnesses and panel bending and twisting stiff- 
nesses. These quantities are not amenable to simple calculation for the complex panels 
tested, and the validity of the values used should be established. 

The effect of K'p on the critical flutter dynamic pressure of panel I is shown in 
figure 12 where A cr is plotted against Kp from 0 to 1.0 and l/Kp from 1.0 to 0. 

The curve represents the critical flutter dynamic pressure over the entire range of Kp. 

In previous investigations of corrugation-stiffened panels, Kp was assumed to be zero 
(for example, ref. 4), since measured natural frequencies were either unavailable or insuf- 
ficient to evaluate the torsional boundary restraint properly. However, as can be seen, 

K'p has a major influence on A cr for panel I and is capable of increasing A cr from 
435 for K'p = 0 to 16 725 for Kip = °°. ^The value of A cr for K-p = 00 is identical to 
the value for Kj) = °°.^ The symbol shown in figure 12 represents the extrapolated exper- 
imental result for panel I and is seen to be only 35 percent of the predicted value based on 
the values of spring constants given in table IV. However, if K-p is neglected (left-hand 
limit of curve), the predicted value of A cr is seen to be about an order of magnitude 
below experiment. It is apparent that if Krp were adjusted from 1.1 to 0.3, theory and 
experiment would be in agreement; however, this value of K'p resulted in theoretical 
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natural frequencies of the panels that differed considerably from those measured. Thus, 
large adjustments in Kq- from the value shown by the symbol are considered to be 
unwarranted. 

A study was also made of the effects of the panel bending stiffness and twisting 
stiffness Dxy on vibration and flutter characteristics of the panels. These stiffnesses 
are of major importance in classical panel flutter analysis but are often difficult to deter- 
mine accurately for a fabricated corrugation-stiffened panel; any inaccuracies in their 
values could cause erroneous flutter results. The results of the study are given in 
table V. 

Table V shows the effect of arbitrary 50 -percent reductions in Dj and D X y on 
flutter dynamic pressure q cr . In addition, changes in several of the panel natural fre- 
quencies are shown as a function of the "best agreement" theoretical natural frequencies. 
(See solid curves in fig. 9.) The reference values of Dj and D X y used in the table 
correspond to the calculated values given in table I. Table V shows that q cr remains 
essentially unchanged for both panels when D X y remains fixed and Dj is reduced. 

Note that the natural frequencies, however, begin to show some disagreement (especially 
for panel I at the higher values of m) with values based on calculated panel stiffnesses 
and spring constants. (The reference frequencies have already been shown in figure 9 to 
be in good agreement with experiment.) Thus, the original calculated value Dj re f is 
adequate for the flutter analysis, since a 50-percent reduction in Dj does not appreciably 
affect q cr . 

In contrast, when Dxy is reduced by a factor of 2 and Dj remains fixed, a corre- 
sponding reduction in q cr of approximately 50 percent is realized. Although this reduc- 
tion in q cr would bring theory and the extrapolated flutter point (for zero in-plane 
loading) into much better agreement, a check of the natural frequencies shows poor agree- 
ment with the reference values for both panels. Even when drastic changes in spring 
constants were explored, the reference frequencies could not be satisfactorily matched. 

It was therefore concluded from this study that it is impossible to predict both the extrap- 
olated flutter points and the natural-vibration characteristics accurately for the test panels 
by using the same set of panel stiffness and support stiffnesses. Inasmuch as more con- 
fidence can be placed in the experimentally determined natural frequencies than the flutter 
points for zero in-plane loading (which are results of large extrapolations), the panel and 
support stiffnesses determined from the vibration analysis are probably reasonably 
accurate. 

Flutter predicti ons for uniform in -plane loading.- Some insight into the discrepancy 
between theory and the extrapolated experimental flutter point may be found by examining 
the theory for panels subjected to uniform in-plane loading. The theoretical flutter 
boundaries for the two panels of this investigation are shown in figures 13(a) and 13(b) 
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in which the cube root of the critical dynamic -pressure parameter is plotted against 
the ratio of uniform in-plane loading in the streamwise direction N x to the buckling 
load N X)Cr . The theoretical boundaries determined from the analysis in appendix B are 
shown by the solid curves for the cross-stream in-plane load Ny = 0. Although the panel 
has a temperature variation through the depth of the corrugations, which sets up a self- 
equilibrating thermal stress system, this stress is neglected in the theory. The exper- 
imental data obtained from figures 6(a) and 6(b) are shown by the symbols. The data are 
plotted by using assumed critical buckling temperature changes near the highest temper- 
ature at which flutter was initiated. The values are 230° F (128 K) for panel I and 80° F 
(44 K) for panel n. Furthermore, it was assumed that in-plane load N x varied directly 
with temperature. Although the experimental data are bracketed by the theoretical values 
of A cr at N x = 0 and Xcr at buckling (N x /N XjCr = l), it is apparent that the data do 
not follow the trend of the theoretical flutter boundaries, especially in the region near 
buckling. It is also apparent that linear extrapolation of the data to zero N x may lead 
to large discrepancies. 

The zero values of the flutter parameter predicted by theory at loads smaller than 
the critical buckling load are characteristic of long, narrow isotropic panels subjected to 
uniform in-plane loads when structural damping is neglected. (See ref. 9.) In fact, on an 
analytical basis, long, narrow isotropic panels are equivalent to short orthotropic panels 
oriented in the airstream as in this investigation. A recent investigation presented in 
reference 10 has indicated that for long, narrow isotropic panels, structural and aero- 
dynamic damping can significantly increase A cr in the region near buckling. This 
increase is illustrated in figure 14, where the cube root of A cr is plotted against the 
uniform in-plane loading ratio for an isotropic panel with length-to-width ratio of 10. 
Experimental flutter data presented in reference 11 for a similar panel are also shown in 
the figure. The theory assumes clamped -edge boundary conditions. The flutter boundary 
obtained with no damping is similar to the boundaries of figure 13 and does not follow the 
experimental trend for values of N x /N x>cr > 0.5. However, with structural and aero- 
dynamic damping included, the trend near buckling is changed dramatically; whereas for 
N x /N X)C r < 0.5, damping has negligible effect on \ cr . Thus, reference 10 indicates that 
when structural damping is included in the theory for isotropic panels, the flutter predic- 
tions near buckling are greatly improved and the theoretical boundary more closely 
follows the experimental trend. 

The influence of structural damping on flutter of spring-supported orthotropic panels 
is beyond the scope of this paper. However, the salient fact is that the amount of experi- 
mental data obtained for the two panels of this investigation is insufficient for linear 
extrapolation of the data as shown in figure 6 to accurately predict flutter dynamic pres- 
sure for N x = 0. In fact, in light of the theoretical boundary for stressed panels, it is 
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believed that structural damping would have to be included in the theory to make a 
meaningful comparison with experiment. 

CONCLUDING REMARKS 

Natural-vibration and flutter characteristics of two corrugation-stiffened panels 
with elastically restrained edges were studied both experimentally and analytically. The 
analysis, which uses small -deflection theory and piston-theory aerodynamics, accounts 
for finite deflectional, rotational, and torsional stiffnesses of the boundary supports and 
uniform in-plane loading. Flutter tests were made for the panels subjected to airflow at 
a Mach number of 3. The panel edges normal to the corrugations (streamwise edges) 
were supported on cantilever beam springs. The support torsional stiffness resulted pri- 
marily from corrugated doubler straps welded to the panel at the ends of the corrugations. 
Panel leading and trailing edges were essentially simply supported. 

The natural -vibration data for the two test panels exhibit excellent agreement with 
theory when the finite deflectional, rotational, and torsional stiffnesses of the boundary 
supports are taken into account in the theory. Only one set of support spring constants 
that permitted correlation between theory and experiment could be found for each test 
specimen. Furthermore, the theory shows that for panels on flexible supports, the addi- 
tion of torsional stiffening at the supports is a definite asset for increasing flutter dynamic 
pressure. This method of support torsional stiffening may be a means for increasing 
flutter resistance of panels for which the supports are designed to permit thermal 
expansion. 

The experimental flutter boundaries obtained as a function of uniform in-plane 
loading exhibit fair agreement with theory at moderate loading levels. Linear extrapo- 
lation of the data to zero uniform in-plane loading resulted in experimental flutter dynamic 
pressures considerably lower than the theoretical predictions. At in-plane loading levels 
near buckling, however, the experimental critical dynamic pressures were much greater 
than theoretical predictions. Examination of theoretical results for the test panels sub- 
jected to uniform in-plane loading indicated that flutter behavior was similar to that of 
panels for which structural damping has been shown to have a pronounced effect at loading 
levels near buckling. It is believed that structural and aerodynamic damping must also be 
included in the theory for flutter of orthotropic panels subjected to uniform in -plane 
loading before a meaningful comparison of theory and experiment can be made. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., August 3, 1970. 
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APPENDIX A 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 

The International System of Units (SI) was adopted by the Eleventh General Confer- 
ence on Weights and Measures, Paris, October 1960. (See ref. 5.) Conversion factors 
for the units used herein are given in the following table: 


Physical quantity 

U.S. Customary 
Unit 

Conversion 

factor 

(*) 

SI Unit 
(**) 

Length 

in. 

0.0254 

meters (m) 

Pressure .... 

lbf/ft 2 

47.88 

newtons/meter 2 (N/m 2 ) 

Stiffness 

lbf-in. 

0.113 

newton-meter (N-m) 

Strength 

lbf/in 2 

6.895 x 10 3 

newtons/meter 2 (N/m 2 ) 

Temperature . . 

°F 

(5/9)(F + 459.67) 

kelvins (K) 


^Multiply value given in U.S. Customary Unit by conversion factor to obtain 


equivalent value in SI Unit. 

**Prefixes to indicate multiples of units are as follows: 


Prefix 

Multiple 

kilo (k) 

10 3 

centi (c) 

10" 2 

milli (m) 

10 -3 
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APPENDIX B 


FLUTTER ANALYSIS OF PANELS ON ELASTIC SUPPORTS 


Flutter Solution 

The orthotropic panel and the coordinate system are shown in figure 7. The leading 
and trailing edges (x = 0 and x = a, respectively) are simply supported; the spanwise 
edges (y = ±j^ are supported by deflectional, rotational, and torsional springs. The 
spring constants per unit length are Kp>, Kr, and Kj corresponding to deflection, 
rotation, and torsion. The panel is loaded by uniform in-plane loads N x and Ny. Air- 
flow at Mach number M passes over one surface of the panel in the x-direction (normal 
to the corrugations). The aerodynamic loading is given by the two-dimensional static 
approximation. 

The present analysis closely parallels the model analysis presented in appendix B 
of reference 3. However, in reference 3 only deflectional springs were considered, and 
the effect of in-plane load was not included. With the addition of the in-plane load terms, 
the governing differential equation given in reference 3 becomes 


Dj ^ + 2 Di 2 

dX* 


dx 2 dy 2 


+ d 2 ^ + N x 

9y’ 8x 2 


+ N t 


8 2 W . o 3 2 W 


+ 2 


2q 


9y^ 


9t 2 Vto2 - 1 


0 


(Bl) 


where 


Di = 


D 2 = 


Dx 

1 ” Mx^y 

Dy 


1 - MxMy 

D12 = D xy + p x D 2 


The boundary conditions at the leading and trailing edges, that is, simple supports, 

are 

w(0,y,r) = w(a,y,r) = -^(0,y,r) = -^(a,y,r) = 0 (B2) 

ax 2 9x 2 

At the spring-supported (streamwise) edges, the boundary conditions to be satisfied can be 
obtained from appendix C of reference 12. The results are 
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APPENDIX B - Continued 





±Kj)W + Krj, 


a 2 w 

9x 2 



9 2 W 
9x 2 9y 


a y 3 


N — 

y ay 


= o 


y=± 2 


(B3) 


The Galerkin procedure is employed to obtain solutions of equation (Bl). A solu- 
tion is assumed to be of the form 


w 


^ A m Y 

m=l X J 


(B4) 


where Y m ^0 is the fundamental cross-stream panel mode shape for no airflow and is 
dependent on m, since the side edges can deflect. By following the procedure of refer- 
ence 3, equation (B4), which satisfies equation (B2), is then substituted into equation (Bl) 
The resulting equation is multiplied by Y n(^) sin(M^) and is integrated over the area 
to obtain the following set of simultaneous equations: 




nn 


J 

^y 


(-i) 


m+n 


2mn 


(n2 - m 2 ) C m,m 


m ’ n A n = 0 


(B5) 


where m = 1, 2, 3, . . j; n / m; and 

1/2 

C m,n - l 1/2 Y m(S Y ngHl) 


(B6) 


The determinant of the coefficients of A m (eq. (B5)) is the flutter determinant which 
must vanish for nontrivial solutions to exist. The frequency w is taken as the eigenvalue. 


Natural Mode Shapes 

The cross-stream natural mode shapes Y m must be such that boundary conditions 
(eqs. (B3)) are satisfied. A suitable set of Y m can be obtained by solving the corre- 
sponding natural vibration problem This solution can be accomplished by substituting 
equation (B4) into equation (Bl) with the aerodynamic loading term set equal to zero. This 
procedure yields the following differential equation for each Y m : 
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APPENDIX B - Continued 


where 


Y"" 

m 


+ 7r 2 Ay jm Y 1 ^ 


^By >m Y m = 0 



(B7) 


(B8) 


and 



d 12 

d 2 


(B9) 


The solution of equation (B7) can be written as 
Y m = Cj^j cosh + Dm cos ^ 


where 


6 


8 










(BIO) 


(Bll) 


Equations (B4) and (BIO) are substituted into the boundary conditions (eqs. (B3)) in 
order to obtain two simultaneous equations in terms of the coefficients C m and D m . 

By setting the determinant of these coefficients equal to zero, a transcendental equation is 
obtained which, when satisfied, insures that the boundary conditions are also satisfied. 

The resulting transcendental equation (for /i x = 0) is 


e 6^ cos tanh |- + 6e ^ sin - irH 
2 > 2 2 


K D + Kt 


(€] 


6 

cos - 
2 


- K R sin | tanh |^(e 2 + 62) 


- 7T 4 K- 


R 


K D + K T /mb ' 


_ . 6 6 , , 
6 sin - + e cos - tanh 
2 Z 


D- 


(B12) 
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APPENDIX B - Concluded 


where 




k» = M 

R £>2ir 


_ K T b 
K T = — =- 
D 2 n 


(B13) 


The expression for the natural mode shape Y m given in equation (BIO) can be rewritten 
as 


V 

Cm 


cosh ^ + — ^2- cos ^ 
b C m b 


(B14) 


where 

D m g2 cos f + ^K r 6 sin | 
Cm e 2 cosh ^ + 7rK R e sinh |- 


(B15) 


Kip on panel 
terms can be thought of as 


An examination of equation (B12) reveals the beneficial effect of large 
flutter dynamic pressure. The sum of the Kp> and j 

an effective spring stiffness term which controls panel flutter dynamic pressure. Thus, 
if the value of Kp can be made large enough, the flutter dynamic pressure of a simply 
supported panel can be approached regardless of how small the value of Kp) may be. It 
should also be noted that Kp becomes more effective with increasing m. 

Equations. (B6), (B8), (B9), (Bll), (B12), (B14), and (B15) are used to determine the 
elements of the flutter determinant in the following manner: 


1. For a given panel, Dp, D 2 , Dp 2 > N x , Ny, and a/b are specified; thus Ay ?m 
(eq. (B9)) is known. 

2. An iteration procedure is used to solve for By jin (eq. (B8)) to determine the 
values of e and 6 which satisfy equations (Bll) and (B12) simultaneously. The solu- 
tion values of e and 6 determine the cross-stream mode shape (eqs. (B14) and (B15)), 
which is used to evaluate C mjn (eq. (B6)). When these steps are completed, the ele- 
ments of the flutter determinant are known, and the flutter behavior of the panel under 
consideration can be determined. 
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TABLE I.- PANEL BENDING AND TWISTING 
STIFFNESS PARAMETERS 

[Calculated by using method of ref. 6] 


Di 7.508 lbf-in. (0.848 N-m) 

D 2 6484 lbf-in. (732.7 N-m) 

Dxy 1342 lbf-in. (151.6 N-m) 


TABLE II.- AVERAGE NATURAL FREQUENCIES 


m 

Panel I 

Panel II 

Average 

frequency, 

Hz 

Deviation, 

Hz 

Average 

frequency, 

Hz 

Deviation, 

Hz 

i 

97 

±2 

127 

±2 

2 

V 

194 

±2 

154 

±2 

3 

290 

±4 

195 

±2 

4 

405 

± 2 

238 

±5 

5 

494 

±4 

294 

±10 

6 

631 

±1 

345 

±4 

7 

752 

±1 

397 

±9 

8 

1 

— 

— 

454 

±5 




TABLE HI.- DATA FOR PANEL FLUTTER TESTS 




Total 


Start of flutter 



End of flutter 




temperature 











Panel 

Test 



^cr 

AT 


Qcr 

AT 




°F 






f, 





f, 



K 












lbf 

kN 

op 

K 

Hz 

lbf 

kN 

Op 

K 

Hz 





ft2 





ft2 

m2 








4931 

236 

133 

74 

205 

4064 

195 

135 

75 

205 


1 

318 

432 

4729 

226 

126 

70 

205 

4307 

206 

110 

61 

205 





4589 

220 

158 

88 

205 

3823 

183 

194 

108 

205 


2 

408 

482 






3432 



109 

205 

I 


4000 

192 

198 

110 

205 

164 

197 





4933 

236 

162 

90 

195 

4094 

196 

213 

118 

200 


3 

452 

506 

3960 

190 

226 

126 

200 







a 4 

360 

455 












1 

308 

426 

2975 

142 

58 

32 

180 







2 

303 

424 

2470 

118 

86 

48 

180 






n 

3 

304 

424 

3230 

155 

36 

20 

255 

2440 

117 

136 

76 

250 


4 

399 

477 

2731 

131 

82 

46 

185 







a No flutter observed. 


TABLE IV.- PANEL SUPPORT SPRING CONSTANTS 







Panel 

k d 

KR 

Kt 

I 

0.918 

0.098 

1.1 

n 

7.1 

2.5 

8 
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TABLE V.- INFLUENCE OF PANEL BENDING STIFFNESS Di AND TWISTING 
STIFFNESS Dxy ON VIBRATION AND FLUTTER BEHAVIOR 


Panel 

D 1 

“ 

D xy 

Qcr 




w m,ref 

(b) 

Dl,ref 

Dxy, ref 


(a) 

(a) 

lbf/in^x 

kN/ m 2 

m = 1 

m = 3 

m = 5 


1 

1 

119.3 

822.6 

1.0 

1.0 

1.0 

I 

i 

1/2 

1 

125.7 

866.7 

1.0 

.97 

.92 

I 

1 

1/2 

68.7 

474 

.90 

.83 

.86 

I 

1 

1 

79.3 

547 

1.0 

1.0 

1.0 

n 

1/2 

1 

80.0 

552 

1.0 

.99 

.96 


1 

1 

1/2 

37.6 

259 

.97 

.80 

.79 


a Di re; £ and D X y re j values are given in table I. 

^Theoretical frequencies w m ,ref based on D-^ re f and D X y, re f. 
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(a) Panel I. 



(b) Panel II. 


Figure 11.- Theoretical flutter mode 
with zero uniform in-plane 


shapes for panels 
loading. 
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